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Compared to healthy controls, individuals with autism spectrum disorders (ASD) have weaker posterior- 
anterior connectivity that strengthens less with age within the default network, a set of brain structures 
connected in the absence of a task and likely involved in social function. The serotonin transporter-linked 
polymorphic region (5-HTTLPR) genotypes that result in lowered serotonin transporter expression are asso- 
ciated with social impairment in ASD. Additionally, in healthy controls, low expressing 5-HTTLPR genotypes 
are associated with weaker default network connectivity. However, in ASD, the effect of 5-HTTLPR on the de- 
fault network is unknown. We hypothesized that 5-HTTLPR's influence on posterior-anterior default network 
connectivity strength as well as on age-related changes in connectivity differs in the ASD group versus con- 
trols. Youth with ASD and healthy controls, ages 8-19, underwent a resting fMRI acquisition. Connectivity 
was calculated by correlating the posterior hub of the default network with all voxels. Triallelic genotype 
was assessed via PCR and Sanger sequencing. A genotype-by-diagnosis interaction significantly predicted 
posterior-anterior connectivity, such that low expressing genotypes (S/S, S/L G , L G /L G ) were associated with 
stronger connectivity than high expressing genotypes (L A /L A , S/L A , L A /L G ) in the ASD group, but the converse 
was true for controls. Also, youth with ASD and low expressing genotypes had greater age-related increases 
in connectivity values compared to those with high expressing genotypes and controls in either genotype 
group. Our findings suggest that the cascade of events from genetic variation to brain function differs in 
ASD. Also, low expressing genotypes may represent a subtype within ASD. 

© 2012 The Authors. Published by Elsevier Inc. All rights reserved. 



1. Introduction 

Autism spectrum disorders (ASD) are neurodevelopmental condi- 
tions characterized by social and communicative impairments and 
rigid repetitive behaviors. The prevalence of ASD has sharply increased 
in recent years and is currently 1 in 88 (CDC, 2012). Deciphering the 
complex etiology of ASD is thus a priority, and progress will likely 
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involve examining the condition using multiple methodologies, includ- 
ing neuroimaging and molecular genetics. 

As alterations in brain connectivity have been repeatedly implicated 
in ASD (Hughes, 2007), attention has been focused on identifying pertur- 
bations in fundamental, large-scale networks, such as the default net- 
work, that may contribute to ASD symptoms. In healthy adults, the 
default network (including the posterior cingulate, angular gyri, superior 
frontal gyri/Brodmann's area (BA) 10, and anterior cingulate/BA 10) 
is active and functionally connected in the absence of a demanding 
task (Raichle and Snyder, 2007). Functional connectivity reflects struc- 
tural connectivity of the default network in healthy adults (Greicius et 
aL, 2009). The default network contains posterior and anterior hubs 
(Buckner et al., 2008) that typically display strong long-range connectiv- 
ity but are distinct from one another (Horovitz et al., 2009). 
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The primary purpose of the default network is a subject of debate. 
The default network may relate to basic central nervous system func- 
tions such as maintaining the balance of excitatory and inhibitory in- 
puts or interpreting information from the environment (Raichle and 
Snyder, 2007). Alternatively, the primary purpose of the default net- 
work may be related to social cognition, including self-referential 
processes (Gusnard et aL, 2001) and mentally projecting oneself 
into hypothetical situations (Buckner and Carroll, 2007). 

Studies on adults with ASD (Cherkassky et al., 2006; Kennedy 
and Courchesne, 2008; Monk et al., 2009) as well as adolescents 
(Weng et al, 2010; Anderson et al., 2011; Wiggins et al., 2011) 
found weaker connectivity between the posterior and anterior de- 
fault network compared to controls. Moreover, the weaker the pos- 
terior-anterior default network connectivity, the worse the social 
impairment in individuals with ASD (Monk et al., 2009; Weng et 
al., 2010). 

A few studies have investigated the development of the default net- 
work. For healthy individuals, posterior-anterior connectivity is weaker 
during childhood and adolescence than during adulthood both func- 
tionally (Fair et al., 2008; Stevens et al., 2009; Wiggins et al., 2011) 
and structurally (Supekar et al., 2010). These studies indicate that con- 
nectivity of this network increases in strength over childhood and ado- 
lescence in healthy individuals. In contrast, youth with ASD have 
attenuated increases in posterior-anterior connectivity with age com- 
pared to controls (Wiggins et al., 2011). 

Identifying the genetic factors that influence the default network in 
ASD is important to further elucidate the complex etiology of ASD. The 
serotonin transporter-linked polymorphic region variant (5-HTTLPR; 
Lesch et al., 1996) in the promoter region of the serotonin transporter 
gene (SLC6A4) is relevant to the default network in ASD. The S and L G 
alleles of 5-HTTLPR are associated with decreased serotonin trans- 
porter expression relative to the L A allele (A to G SNP in the L allele, 
rs25531; Hu et al., 2006). The low expressing alleles of 5-HTTLPR 
have been associated with worse social symptoms in ASD (Tordjman 
et al., 2001; Brune et al., 2006). In healthy adolescents, 5-HTTLPR is 
known to influence the default network: those with low expressing 
genotypes exhibit weaker posterior-anterior connectivity than ado- 
lescents with high expressing genotypes (Wiggins et al., 2012). More- 
over, in healthy children and adolescents, 5-HTTLPR also impacts the 
development of default network connectivity such that youth with 
high expressing genotypes have greater age-related increases in pos- 
terior-anterior connectivity than those with low expressing geno- 
types (Wiggins et al., 2012). A previous study found that serotonin 
transporter binding in the anterior default network is decreased in in- 
dividuals with autism (Nakamura et al., 2010). However, no study has 
yet examined how 5-HTTLPR affects default network connectivity or 
its development in individuals with ASD. 

The present study addresses these two gaps in the literature on ASD: 
the role of 5-HTTLPR in default network connectivity and in the develop- 
ment of default network connectivity. This is accomplished by directly 
examining the influence of 5-HTTLPR variants on posterior-anterior de- 
fault network connectivity as well on as age-related changes in connec- 
tivity in a sample of children and adolescents with ASD and controls. We 
hypothesized that the relationship between 5-HTTLPR genotype and 
posterior-anterior default network connectivity strength differs in the 
ASD group versus controls. Additionally, we hypothesized that the rela- 
tionship between 5-HTTLPR and changes in connectivity across child- 
hood and adolescence differs in the ASD group compared to controls. 

2. Material and methods 

2.1. Participants 

Fifty-four children and adolescents with ASD and 66 healthy con- 
trols, aged 8.3 to 19.6 years, were included in this study (see Table 1 
for participant characteristics). From a total of 105 participants with 



ASD and 82 controls recruited, 51 participants with ASD and 16 con- 
trols were excluded because of head movement exceeding 2.5 mm 
translation or 2.5° rotation, declining to complete the MRI scan due 
to discomfort, failure to return a saliva sample for genotyping, or 
technical problems with the MRI. 

Controls were recruited through flyers posted at community organi- 
zations in the Ann Arbor, Michigan area. The University of Michigan 
Autism and Communication Disorders Center (UMACC) referred poten- 
tial participants to our study and diagnosed participants with an ASD 
(autistic disorder, Asperger's syndrome, or pervasive developmental dis- 
order — not otherwise specified) using the Autism Diagnostic Obser- 
vation Schedule (ADOS; Lord et al., 2000), the Autism Diagnostic 
Interview— Revised (ADI-R; Lord et al., 1994), and clinical consensus 
(Lord et al., 2006). The University of Michigan Institutional Review 
Board approved the procedures. Participants over age 18 gave written 
informed consent; participants under age 1 8 gave written assent and 
their parents gave written informed consent. Cognitive functioning 
was evaluated for controls with the Peabody Picture Vocabulary Test 
(PPVT; Dunn and Dunn, 1997) and the Ravens Progressive Matrices 
(Raven, 1960); participants with ASD were given these measures or 
the Differential Ability Scales II - School Age (Elliott, 2005), the 
Stanford-Binet Intelligence Scales (Roid, 2003), the Wechsler Intelli- 
gence Scale for Children IV (Wechsler, 2003), or the Wechsler Abbre- 
viated Scale of Intelligence (Wechsler, 1999). Participants with 
orthodontic braces, medical conditions contraindicated for MRI, or his- 
tory of seizures or neurological disorders were excluded. Control partic- 
ipants were screened for psychological disorders with the Child 
Behavior Checklist (Achenbach and Edelbrock, 1981), Social Respon- 
siveness Scale (Constantino et al., 2003), Social Communication Ques- 
tionnaire (Rutter et al., 2003), Obsessive Compulsive Inventory — 
Revised (Foa et al., 2010), Child Depression Inventory (Kovacs, 1992), 
and Multidimensional Anxiety Scale for Children (March et al., 1997). 
All control participants scored below clinical cutoffs for affected status. 
Individuals with the low and high expressing genotypes did not differ 
in any of the symptom measures or cognitive functioning in either the 
ASD or control group ( Inline Supplementary Table SI ). Prior studies uti- 
lized portions of this dataset (Weng et al., 2010, 2011; Wiggins et al., 
2011,2012). 

Inline Supplementary Table SI can be found online at http://dx. 
doi.org/10.1016/j.nicl.2012.10.008. 



2.2. Genetic analyses 

5-HTTLPR genotype was ascertained using previously published 
procedures (Wiggins et al., 2012). Participants donated saliva samples 
using the Oragene DNA kit (DNA Genotek; Kanata, Canada). PCR and 
agarose genotyping were used to determine S versus L allele. Sanger se- 
quencing was utilized to determine the A to G single nucleotide poly- 
morphism (SNP) in the L allele (rs25531; Hu et al., 2006) and to 
confirm PCR genotyping. 

In autism, individuals with the low expressing genotype (S/S) have 
been shown to differ in neurochemical metabolism compared to L allele 
carriers in the anterior portion of the default network (Endo et al., 
2010). As such, participants were put into two genotype groups: low 
expressing genotypes (S/S, S/Lq, Lq/Lg) versus medium and high ex- 
pressing genotypes (La/L a , S/L a , La/Lg, hereafter referred to as "high ex- 
pressing" genotypes). (The Lg allele is equivalent to the S allele in 
serotonin transporter expression level (Hu et al., 2006), so for the pur- 
poses of the analyses, the two alleles were grouped together.) This ge- 
notype grouping is consistent with a number of non-ASD studies that 
found recessive effects of the low expressing 5-HTTLPR alleles, often 
in adolescent populations (e.g., Cicchetti et al., 2007; Surguladze et al., 
2008; Benjet et al., 2010). Nevertheless, we conducted additional anal- 
yses to examine whether our results still stood when the alleles were 
grouped differently (see Section 3.1.4 Alternative genotype groupings). 
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Hardy-Weinberg equilibrium was tested based on the insertion/ 
deletion polymorphism. Genotype frequencies were in Hardy-Weinberg 
equilibrium for the ASD group {/ 2 = 0.742, df= 1, p = 0.389), but there 
was a trend toward disequilibrium for the control group (j 2 = 3.74, df = 
1 , p = 0.053). When including only Caucasians for the control group, the 
trend disappeared {x 2 = 0.654, df= 1, p = 0.419). Because of this, addi- 
tional analyses were performed to address the potential effects of multi- 
ple ancestries within the sample. 

2.3. JMRI data acquisition 

T 2 *-weighted blood oxygen level dependent (BOLD) images were ac- 
quired during a resting state scan, in which participants were instructed 
not to think of anything in particular and to let their minds wander while 
looking at a fixation cross. Over the 10-minute resting state scan, 300 im- 
ages were acquired (Glover and Law, 2001; TR=2000 ms, TE = 30 ms, 
flip angle = 90°, FOV = 22 cm, 64 x 64 matrix, 40 contiguous axial 
3 mm slices). Slices were acquired parallel to the intercommissural 
(AC-PC) line. High-resolution 3D Tl axial overlay (TR = 8.9, TE=1.8, 
flip angle = 15°, FOV = 26 cm, slice thickness = 1.4 mm, 124 slices; ma- 
trix=256xl60) and spoiled gradient (SPGR; flip angle = 15°, FOV= 
26 cm, 1.4 mm slice thickness, 110 slices) anatomical images were also 
collected. Participants wore a pulse oximeter and abdominal pressure 
belt to record cardiac and respiratory rhythms, synchronized to the 
fMRl data, for subsequent physiological artifact correction. Further details 
on the acquisition parameters have been previously published (Weng et 
al., 2011; Wiggins et al., 2012). Prior to the MRI scan, participants prac- 
ticed in a mock scanner to acclimate to the scanning environment. 

2.4. fMRI data analysis 

2.4.1. Data preprocessing 

The standard pre-processing procedure from the University of Mich- 
igan Functional MRI Center was applied to the fMRI data. This procedure 
includes removing outliers from the raw k-space data, reconstructing 
the k-space data to image space, applying a field map correction to re- 
duce artifacts from susceptibility regions, and correcting for slice timing. 
RETROICOR was utilized to remove noise associated with cardiac and 
respiratory rhythms (Glover et al., 2000). To address potential effects 
of head motion, functional images were realigned to the 10th image. 
Details on these steps are available in multiple papers utilizing this 
pre-processing stream (e.g., Weng et al., 2011; Wiggins et al., 2012). 
The high-resolution Tl anatomical images were then co-registered to 
the functional images using the SPM5 Matlab toolbox (Wellcome 
Department of Neurology, London, UK; http://www.fil.ion.ucl.ac.uk). 
After removing variance associated with head motion (see Section 
2.4.2 Addressing head motion), functional images were smoothed 
with an isotropic 8 mm full width at half maximum (FWHM) Gaussian 
kernel using SPM5. A low-pass filter of .08 Hz was applied as well to iso- 
late the frequency band where default network activation has been 
found. 

2.4.2. Addressing head motion 

Recent papers have emphasized the importance of addressing 
head motion, which can introduce spurious correlations in connectiv- 
ity analyses (Power et al., 2012; Satterthwaite et al., 2012; Van Dijk et 
al., 2012). In addition to the standard realignment of images to cor- 
rect for head motion, we took several steps to address the potential 
effects of head motion on our results. First, we excluded participants 
whose head motion exceeded 2.5 mm in the x, y, or z direction or 
2.5° in the roll, pitch, or yaw directions. 

Second, we removed variance associated with head motion by cre- 
ating nuisance regressors from motion estimated in the x, y, z, roll, 
pitch, and yaw directions and retaining the residuals for processing. 

Third, we calculated the mean motion for each person (i.e., mean ab- 
solute displacement of each volume as compared to the previous volume, 



calculated as square root ( (x,+ ! — x f ) 2 + (y,+ 1 — yO 2 + (z, +1 — z,) 2 ) for 

i = l 300 images), as in Van Dijk et al (2012). We then conducted a 

2-way ANOVA with post-hoc contrasts to examine whether mean mo- 
tion differed between the ASD and control groups and by genotype. We 
also performed a 3-way interaction analysis to examine whether the in- 
teraction of genotype-by-diagnosis-by-age significantly predicted mean 
motion. 

Fourth, even though groups may not differ on overall mean motion, it 
is possible that the different motion distributions within groups could in- 
fluence findings. Because of this, we repeated the analyses for our main 
hypotheses with a subsample of our participants matched on mean mo- 
tion to examine whether the results persisted when the groups' motion 
distributions were equivalent (see Section 3.1.1 Matched head motion 
distributions). Several other studies have utilized matching to reduce 
the likelihood that the connectivity patterns they observed were an arti- 
fact of head motion (Fair et al., 2007, 2008; Dosenbach et al., 2010). Our 
matching procedure was as follows: first, we split participants into four 
groups: individuals with the low expressing genotypes (ASD and con- 
trol) and individuals with the high expressing genotypes (ASD and con- 
trol). Within each group, we binned participants by mean motion into 
.001 mm bins. Participants were removed randomly until the number 
of participants in each corresponding bin for the ASD and control groups 
were the same within both the low and high expressing genotypes. 

2.4.3. Connectivity images 

A self-organizing map algorithm was applied to the images to derive 
a data-driven reference from the posterior hub of the default network to 
calculate connectivity for each individual, as described in previous pub- 
lications (Peltier et al., 2003; Wiggins et al., 2011, 2012). An example of 
a posterior hub identified for an individual using the self-organizing 
map algorithm is shown in Inline Supplementary Fig. SI . The advantage 
of using this data-driven method is that, unlike traditional a priori seed 
analyses, seed placement is not based on data from adult control brains. 
When using the self-organizing map algorithm, the default network ref- 
erence, which is correlated with every other voxel in the brain to calcu- 
late connectivity, is not biased toward the control group but rather 
tailored for each individual. 

Inline Supplementary Fig. SI can be found online at http://dx.doi. 
org/10.1016/j.nicl.2012.10.008. 

The connectivity images generated by this method were normalized 
to Montreal Neurological Image ( MNI) space by estimating the transfor- 
mation matrix for the SPGR image to SPM's template MNI image, and 
applying that transformation to the connectivity images. The end prod- 
uct is a normalized image for each subject that indicates, with a Z value 
at each voxel, how highly functionally connected ( correlated ) that voxel 
is to the posterior hub of the default network identified by the 
self-organizing map algorithm. 

2.4.4. Group-level analyses 

The connectivity images were then entered into second-level analy- 
ses in SPM8 to test hypotheses at a group level. As a preliminary step, 
we first examined whether the ASD group and the control group 
exhibited default network connectivity by applying small volume cor- 
rections using masks covering the default network (posterior cingulate, 
precuneus, angular gyri, inferior parietal lobules, parahippocampal gyri, 
superior frontal gyri, anterior cingulate, BA 32, BA 10), from the Wake 
Forest Pickatlas (Maldjian et al., 2002). Also as a preliminary step, we 
compared the ASD and control groups on long-range default network 
connectivity, using the anterior masks in small volume corrections 
(Inline Supplementary Table S2). 

Inline Supplementary Table S2 can be found online at http://dx. 
doi.org/1 0.1 01 6/j.nicl.201 2.1 0.008. 

To address our first hypothesis, a voxel-wise multiple regression 
was created to examine the interaction of genotype (low expressing 
(S/S, S/Lc, Lg/Lg) versus high expressing (La/L a , S/L a , L a /Lg)) by diagno- 
sis (ASD versus control group). For this model, three regressors were 
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entered - genotype, diagnosis, and the interaction of genotype-by- 
diagnosis - predicting connectivity with the posterior hub. To deter- 
mine whether the beta for the interaction was significant in the anterior 
default network, a small volume correction was performed on the 
image mapping the betas of the interaction using a mask of the bilateral 
BA 10, where alterations in long-range default network connectivity 
have consistently been found in ASD samples (e.g., Monk et al., 2009; 
Weng et al., 2010; Wiggins et al., 2011), as well as where effects of 
5-HTTLPR have been found (Wiggins et al., 2012). The small volume 
correction takes into account the geometric qualities of the mask 
when doing a correction for multiple comparisons based on the number 
of resels (a measure related to the number of independent observa- 
tions) within the mask (Worsley et al., 1996). Significance thresholds 
within BA 10 were corrected for multiple comparisons using family 
wise error (FWE) correction (Worsley et al., 1996). Post-hoc compari- 
sons were also performed in SPSS, comparing each subgroup pair 
on connectivity values extracted and averaged from a 4 mm sphere 
around the peak of the interaction, with a Bonferroni-corrected a 
level of 0.05/6 = 0.0083. 

To address our second hypothesis, we created a model to examine 
the three-way interaction among genotype (low versus high express- 
ing), diagnosis (ASD versus control group), and age. In this model, the 
three-way interaction term was entered, as well as all lower order 
terms. A small volume correction applied the same mask as in the first 
hypothesis, bilateral BA 10, to the image mapping the betas for the 
3-way interaction to examine whether the three-way interaction signif- 
icantly predicted connectivity with the posterior hub in the anterior de- 
fault network. Post-hoc analyses were also performed to further 
characterize the interaction. Connectivity values from a 4 mm sphere 
around the peak of the 3-way interaction were extracted and averaged, 
then exported to SPSS. In SPSS, the simple slopes (changes in connectiv- 
ity for every unit increase in age) for four subgroups (controls with low 
expressing genotypes, controls with high expressing genotypes, indi- 
viduals with ASD and low expressing genotypes, individuals with ASD 
and high expressing genotypes) were tested against zero. 

3. Results 

Individuals with the low and high expressing genotypes within 
the control and the ASD groups did not differ on any of the symptom 
measures (Inline Supplementary Table SI). Both the ASD group and 
the control group exhibited default network connectivity, and previ- 
ous findings of weaker posterior-anterior default network connectiv- 
ity in the ASD group were replicated ( Inline Supplementary Table S2; 
Cherkassky et al., 2006; Kennedy and Courchesne, 2008; Monk et al., 
2009; Weng et al., 2010; Wiggins et al., 2011). 

The four groups (individuals with ASD and low expressing geno- 
types, individuals with ASD and high expressing genotypes, controls 
with low expressing genotypes, controls with high expressing geno- 
types) did not differ in mean head motion (genotype-by-diagnosis: 
Fi,ii6 = 0.040, p = .841 ). Additionally, age did not relate to head motion 
differently across the four groups (genotype-by-diagnosis-by-age; p = 
.058, tn 2 = 0.256, p = 0.799). 

The first hypothesis, that the relationship between 5-HTTLPR geno- 
type and posterior-anterior default network connectivity strength dif- 
fers in the ASD group versus controls, was confirmed. There was a 
significant genotype-by-diagnosis interaction predicting degree of con- 
nectivity between the posterior hub and the anterior default network 
in the left hemisphere (xyz=-34, 62, 0, r 116 = 4.24, p = 0.021, 
corrected for multiple comparisons within bilateral BA 10; Fig. 1). Spe- 
cifically, 5-HTTLPR genotype influences posterior-anterior connectivity 
strength differently for individuals with ASD versus controls. Two 
pair-wise comparisons survived a Bonferroni correction, indicating 
that individuals with low expressing genotypes within the ASD group 
had significantly stronger connectivity than individuals with ASD and 
high expressing genotypes (p = 0.001 ) as well as controls with low 



expressing genotypes (p = 0.003). The genotype-by-diagnosis interac- 
tion was also significant in the right anterior default network (xyz = 
44, 56, —6, r 116 = 4.17, p = 0.027, corrected for multiple comparisons 
within bilateral BA 10; Inline Supplementary Fig. S2). 

Inline Supplementary Fig. S2 can be found online at http://dx.doi. 
org/1 0.1016/j.nicl.2012.10.008. 

Our second hypothesis, that the relationship between 5-HTTLPR and 
changes in posterior-anterior connectivity across childhood and adoles- 
cence differs in ASD compared to controls, was also confirmed. We 
found a significant genotype-by-diagnosis-by-age interaction predicting 
degree of connectivity between the posterior and anterior default net- 
work (xyz = — 6, 40, — 6, r 112 = 4.09, p = 0.037, corrected for multiple 
comparisons within bilateral BA 10; Fig. 2). 5-HTTLPR genotype differ- 
entially influences age-related changes in posterior-anterior con- 
nectivity strength in individuals with ASD compared to controls. 
Post-hoc analyses to further characterize the interaction indicated that 
only individuals with ASD with low expressing genotypes had signifi- 
cant increases in connectivity values with age (simple slope = 0.708, 
p = 0.002), whereas the other subgroups' relationships between con- 
nectivity and age did not significantly differ from zero: ASD group, 
high expressing genotype, simple slope = — .268, p = 0.104; controls, 




Controls ASD 



Fig. 1. Impact of 5-HTTLPR genotypes on posterior-anterior default network connectivity 
is different in youth with ASD compared to controls. Voxels in color indicate places where 
connectivity between that area and the posterior default network is differentially 
influenced by 5-HTTLPR in the ASD group versus controls. A significant genotype-by- 
diagnosis interaction in the anterior default network (xyz= —34, 62, 0, t 116 = 4.24, p = 
0.021 , corrected for multiple comparisons within bilateral BA 1 0) is depicted in the trans- 
verse section of the brain (upper). For this and the subsequent brain image, the threshold 
was set atp<0.01 for illustration purposes. To show the interaction, contrast values from a 
4 mm sphere around the peak voxel (xyz = — 34, 62, 0) were extracted and plotted 
(lower). In the bar graph, controls show the pattern found in previous research (Wiggins 
et al., 201 2), but youth with ASD show a different pattern. Brackets with asterisks indicate 
significant differences at a Bonferroni-corrected a-level of 0.05/6 = 0.0083. 
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low expressing, simple slope = 0.21 6, p = 0.335; controls, high express- 
ing, simple slope = — .185, p = 0.229. 

3.1. Additional analyses 

In imaging and genetic studies with disordered populations, head 
motion, population stratification, psychotropic medication status, allele 
grouping, and degree of smoothing are potential factors influencing as- 
sociations. Because of this, additional analyses were performed to deter- 
mine whether these factors account for our results. As these additional 
analyses required a reduced sample size and/or more complex models, 
thereby diminishing the power to detect effects, we utilized a threshold 
ofp<0.05 without family- wise error correction. To summarize, the hy- 
potheses were still confirmed even when taking into account each of 
these factors. 

3.3.1. Matched head motion distributions 

Our matching procedure is described in Section 2.4.2 Addressing 
head motion. In total, 24 participants (20%) were removed in order 
to match the groups' motion distributions. (See Inline Supplementary 
Fig. S3 for a visual representation of the participants removed from 
each bin.) After removing participants until the group distributions 
were matched on mean head motion, our first hypothesis was still 
confirmed. The interaction of genotype-by-diagnosis significantly 
predicted connectivity in bilateral BA 10 (left: xyz= — 34, 62, 0, 



t 92 = 4.22, p = 0.000029; right: xyz = 38, 60, 2, t 92 = 3.84, p = 
0.0001 1 ) with a subsample of participants matched on head motion. 
Our second hypothesis was confirmed as well with the matched sub- 
sample, as a genotype-by-diagnosis-by-age interaction significantly 
predicted connectivity in BA 10 (xyz= — 10, 36, —6, t 88 = 4.22, p = 
0.000029). 

Inline Supplementary Fig. S3 can be found online at http://dx.doi. 
org/10.1016/j.nicl.2012.10.008. 

3.1.2. Population stratification 

To determine whether the findings were primarily driven by differ- 
ing ancestries within our sample, 5 non-Caucasian individuals with ASD 
and sixteen non-Caucasian controls were excluded and the group-level 
analyses addressing our hypotheses were repeated. In line with our first 
hypothesis, including Caucasian participants only, the genotype-by- 
diagnosis interaction predicting connectivity strength was significant 
in both left (xyz=-34, 62, -2, t 95 = 3.95, p = 0.000075) and right 
(xyz = 38, 58, 4, t 95 = 3.14, p = 0.0011) BA 10. Supporting our second 
hypothesis, the genotype-by-diagnosis-by-age interaction predicting 
connectivity was significant in BA 10 when including only Caucasian 
participants (xyz = - 20, 50, 2, t 91 = 3.51, p = 0.00035). 

3.J.3. Medication effects 

Next, twenty-four individuals with ASD taking psychotropic medi- 
cation and 1 control taking levothyroxine (a thyroid medication) were 
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Fig. 2. 5-HTTLPR influences age-related changes in posterior-anterior default network connectivity differently in youth with ASD compared to controls. Voxels in color indicate 
places where connectivity between that area and the posterior hub changes across age differently for the ASD group and the control group. A significant 
genotype-by-diagnosis-by-age interaction in the anterior default network (xyz= — 6, 40, —6, tn 2 = 4.09, p = 0.037. corrected for multiple comparisons within bilateral BA 10) 
is depicted in the transverse section of the brain (upper). To illustrate connectivity levels in each individual, contrast values from a 4 mm sphere around the peak voxel 
(xyz = — 6, 40, —6) were extracted and plotted (lower). 
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excluded before repeating the analyses. Like the analyses only including 
Caucasians, findings including only non-medicated participants also 
mirrored the original findings with the entire dataset. Supporting the 
first hypothesis, with only non-medicated participants, the genotype- 
by-diagnosis interaction was significant ( left BA 1 0: xyz = — 34, 62, 0, 
t 91 =3.71, p = 0.00018; right BA 10: xyz = 32, 60, 6, t 91 =4.55, p = 
0.0000084). The second hypothesis, genotype-by-diagnosis-by-age in- 
teraction predicting connectivity, also held with non-medicated partic- 
ipants only in BA 10 (xyz = - 6, 40, - 6, t 87 = 2.97, p = 0.0019). 

3. 1.4. Alternative genotype groupings 

Moreover, additional analyses were conducted with alternate geno- 
type groupings to examine whether the findings persisted when partic- 
ipants split into different genotype groups. First, we ran the analyses 
with the genotype in three groups based on expressing level: low ex- 
pressing (S/S, S/Lg, Lg/Lc) versus medium expressing (S/L A , La/Lq) vs 
high expressing (La/La) genotypes. Consistent with the first hypothesis, 
there was a significant genotype-by-diagnosis interaction in both the 
left (xyz= -34, 62, 0, F 2 ,ii4 = 9.22, p = 0.00019) and right (xyz = 44, 
56, -6, F 2114 = 9.19, p = 0.00020) anterior default network. The 
three-way genotype-by-diagnosis-by-age interaction was also signifi- 
cant with this genotype grouping in the anterior default network 
(xyz =-6, 40, -8, F 2108 = 11.23, p = 0.000037), consistent with the 
second hypothesis. 

Second, we examined the hypotheses with the genotype grouping S/S 
versus heterozygotes (S/L A and S/Lc) versus La/L a . With this alternative 
genotype grouping, there was a significant genotype-by-diagnosis in- 
teraction in the left (xyz=-34, 62, -2, F 2114 = 7.68, p = 0.00074) 
and right (xyz = 44, 56, -6, F 2 ,n 4 = 10.45, p = 0.000068) anterior de- 
fault network. Moreover, consistent with the second hypothesis, the 
three-way interaction was significant in the anterior default network 
(xyz=-6, 40, -8, F 2108 = 8.11, p = 0.00052) with this alternative 
grouping. To summarize, the original result pattern persisted even 
when genotypes were grouped two alternate ways in the statistical 
analyses. 

3.1.5. 5 mm smoothing kernel 

The degree of smoothing can also affect results. We re-did the anal- 
yses with a 5 mm (instead of 8 mm) FWHM Gaussian kernel for spatial 
smoothing of the functional images. Consistent with the first hypothe- 
sis, the genotype-by-diagnosis interaction was significant in the left 
(xyz = -48, 54, 6, t 116 = 2.25, p = 0.013) and right (xyz = 42, 58, 10, 
t n6 = 3.96, p = 0.00065) BA 10. Moreover, consistent with the second 
hypothesis, the three-way interaction was significant in BA 10 (xyz = 
42, 48, 24, tn 2 = 3.03, p = 0.002) with the 5 mm smoothing kernel, al- 
beit in a more lateral location within BA 10. 

4. Discussion 

In results confirming both the first hypothesis (genotype-by- 
diagnosis) and the second hypothesis (genotype-by-diagnosis-by-age), 
individuals with ASD and low expressing genotypes stood out among 
the other subgroups, exhibiting the greatest connectivity as well as 
the sharpest increase in connectivity values with age. This overall pat- 
tern suggests that individuals with low expressing genotypes may rep- 
resent a subtype of ASD, which is consistent with previous research 
linking 5-HTTLPR to symptom subtypes rather than a global ASD diag- 
nosis in a larger sample (Tordjman et al., 2001; Brune et al., 2006). 
Linking genotype to brain phenotypes may be a more sensitive way to 
identify subtypes of ASD than linking genotype to behavior, as individ- 
uals with ASD and the low and high expressing genotypes did not differ 
on any of the symptom measures (Inline Supplementary Table S2). Fu- 
ture research could examine other aspects of this potential subtype of 
ASD, including responsiveness to specific interventions and long-term 
prognosis. 



There are two main possibilities to explain why 5-HTTLPR influences 
the ASD group and control group differently. First, a gene-by-gene inter- 
action may account for our results. The specifics of the complex genetic 
etiology of ASD are a subject of intense inquiry. Nonetheless, as ASD is 
highly heritable (Miles, 2011), individuals with ASD may cany a sys- 
tematically different overall genetic profile than controls. Causative 
gene products may interact with 5-HTTLPR, leading to alterations in ex- 
pression levels that then produce a different brain phenotype. Future 
research probing this possibility may include examining other autism 
genes and their involvement in serotonin metabolism. 

Alternatively, a gene-by-environment interaction may explain our 
findings. As 5-HTTLPR is sensitive to environmental input (Belsky et 
al., 2009), it may be that 5-HTTLPR affects brain function in individuals 
with ASD differently than controls because individuals with ASD expe- 
rience an altered social environment brought about by the reactions of 
others to their symptoms. Particularly during adolescence, an important 
social development period in which relationships with peers become 
more important (Youniss and Haynie, 1992), individuals with ASD 
may miss out on social opportunities with peers and thus find them- 
selves in an environment with reduced social stimuli. This environment 
could affect epigenetically-sensitive serotonin transporter expression 
and subsequently, brain function. Future studies incorporating compre- 
hensive environmental measures and focusing on molecular mecha- 
nisms of altered serotonin transporter expression, such as methylation, 
will be necessary to probe this possibility. 

This study has several limitations, as some confounds make imaging 
and genetics research, particularly with pediatric clinical populations, 
more challenging. First, in our study, mean head motion did not differ 
among the groups we compared: individuals with ASD and low ex- 
pressing genotypes, individuals with ASD and high expressing geno- 
types, controls with low expressing genotypes, and controls with high 
expressing genotypes. Also, age did not relate to head motion different- 
ly across the four groups. Nevertheless, motion remains a concern in all 
functional connectivity studies, so we took several steps to address 
motion: first, only participants with movement under 2.5 mm and 
2.5° in all translation and rotation directions were included; second, 
we realigned the functional images; third, we removed variance associ- 
ated with movement in the x, y, z, roll, pitch, and yaw directions; fourth, 
we repeated the analyses with a subsample matched on head motion 
and found that our hypotheses were still confirmed even when motion 
distributions were the same across individuals with ASD and controls in 
the low and high expressing genotype groups. 

Another limitation is that the cross-sectional design utilized in this 
study precludes inferences about developmental trajectories within in- 
dividuals. Future studies may use a longitudinal design to rule out birth 
cohort effects. Additionally, longitudinal studies will be useful to exam- 
ine whether brain differences earlier in development predict later 
symptom presentation and responsiveness to particular treatments. 

Third, we did not exclude any racial or ethnic group when recruiting 
participants, which can contribute to spurious associations in genetic 
studies. Although it should be acknowledged that Caucasians in our 
sample may not all be of the same ancestry, we repeated the analyses 
with non-Caucasian participants removed to determine whether results 
were primarily due to several different ancestries within the sample. The 
genotype-by-diagnosis and genotype-by-diagnosis-by-age interactions 
predicting posterior-anterior connectivity were found even with non- 
Caucasians excluded from the analyses, suggesting that the results 
were not primarily driven by population stratification. Nevertheless, 
the lack of understanding of genetic effects in different racial/ethnic 
groups is a widespread problem in the field that must be addressed in 
future work. 

4.1. Conclusions 

This is the first study, to our knowledge, to examine the influence 
of 5-HTTLPR genotype on the default network in individuals with 
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ASD. We found that the relationship between 5-HTTLPR genotype and 
posterior-anterior default network connectivity is different in indi- 
viduals with ASD compared to controls. Specifically, consistent with 
previous research in controls (Wiggins et aL, 2012), high expressing 
genotypes were associated with stronger connectivity than low ex- 
pressing genotypes. However, the pattern was reversed for the ASD 
group: individuals with ASD and low expressing genotypes had stron- 
ger connectivity than individuals with ASD and high expressing ge- 
notypes. Also, we found that youth with ASD and low expressing 
genotypes had greater age-related increases in connectivity values 
compared to others in the ASD group with high expressing genotypes 
and to controls with either low or high expressing genotypes. The 
present findings provide evidence that the cascade of events from ge- 
netic variation to brain function is markedly different in ASD versus 
typically developing, healthy individuals. Moreover, the findings sug- 
gest that the impact of genotype on brain function is not static but 
rather develops and changes with age. Thus, understanding how 
5-HTTLPR affects brain function in ASD is dependent on the develop- 
mental timeframe. 

Although replication of our findings with a larger sample is neces- 
sary, the present study lays the groundwork to better understand the 
genetic and brain mechanisms that are involved in ASD. The present 
study documented a different impact of 5-HTTLPR on both default net- 
work connectivity and the development of default network connectivi- 
ty in ASD compared to controls. Future studies may expand on these 
findings by examining the structural connections within the default 
network in vivo using diffusion tensor imaging. Moreover, the resting 
connectivity approach used in this study will be useful to examine the 
brain activation patterns of lower functioning individuals with ASD or 
very young children. These individuals are underrepresented in func- 
tional MR1 studies because they are often unable to comply with the de- 
mands of a task requiring responses in the scanner. The relatively low 
demand of a resting fMRl acquisition, on the other hand, may allow 
lower functioning and younger participants to be successfully scanned. 
Obtaining brain data from individuals with a greater range of cognitive 
abilities and ages will allow researchers to gain a broader, more repre- 
sentative picture of ASD and the developmental trajectory of ASD earlier 
than mid-childhood. To conclude, the findings from our study open a 
path for a research program to better understand genetic influences 
on brain function in ASD. 
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